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This work explores the use of functionally graded materials for the aeroelastic tailoring 
of a metallic cantilevered plate-like wing. Pareto trade-off curves between dynamic stability 
(flutter) and static aeroelastic stresses are obtained for a variety of grading strategies. A key 
comparison is between the effectiveness of material grading, geometric grading (i.e., plate 
thickness variations), and using both simultaneously. The introduction of material grading 
does, in some cases, improve the aeroelastic performance. This improvement, and the 
physical mechanism upon which it is based, depends on numerous factors: the two sets of 
metallic material parameters used for grading, the sweep of the plate, the aspect ratio of the 
plate, and whether the material is graded continuously or discretely. 
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Nomenclature 

aerodynamic influence coefficients 
reference length 

force vector due to angle of attack 

failure metric 

damping 

imaginary number 
reduced frequency 
stiffness matrix 

continuous and discrete material fractions 

wing mass and target wing mass 

mass matrix 

aeroelastic eigenvalue 

dynamic pressure 

static solution vector 

critical speed 

flutter and divergence speeds 
aeroelastic eigenvector 
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V = flow speed 

p = flow density 

a) = angular frequency 

[O] = modal matrix 

<r VM , °y = von Mises and yield stress 

I. Introduction 

A potential enabling technology for the aeroelastic tailoring of wing structures may be functionally graded 
materials (FGM), which have received relatively little attention. Functionally graded materials have 
continuously varying properties by spatially varying the distribution of two (or more) materials, which facilitates 
designs with tailored stiffness within a continuous metallic structure [1]. New manufacturing processes, such as 
electron beam freeform fabrication (EBF 3 ), an additive manufacturing process [2], are helping to enable the 
fabrication of functionally graded metals and making FGMs accessible to aircraft design. FGMs offer two 
potentially advantageous design capabilities. First, they enable continuous changes in material properties (elastic 
modulus, density, yield stress, etc.) throughout a structure, allowing local properties to be tuned. An example of this 
is shown in Figure 1, where a stiffener is composed of a blend of stainless steel and copper, and grading can be 
accomplished along the stiffener axis, or through the depth. Second, they enable changes in structural stiffness 
without necessarily requiring a geometric change in the structure, such as an increase in thickness. 



Figure 1. Cross sections of two stiffeners made from a 50/50 SS-Cu blend. Left: grading along the stiffener 
(axis coming out of the page): cut taken mid-length of the stiffener. Right: grading through the depth (axis is 

vertical). Image courtesy of K. Taminger, NASA LaRC. 

In spite of their potential to improve the state-of-the-art in high-performance aircraft structures, few examples of 
FGMs in an aerospace setting are available in the literature [3]. Most of the existing work details the use of FGM 
for elastic panels subject to supersonic flows (and thus, aerodynamic heating): see Refs. [4], [5], and [6]. Gradual 
material grading is particularly beneficial to high-temperature applications since it can eliminate discrete changes in 
the coefficient of thermal expansion, which would introduce stress concentrations at material boundaries [6]. FGMs 
have been shown in these papers to improve the aerothermoelastic panel flutter boundaries as well. For subsonic 
structures, Librescu and Maalawi [7] use material grading to optimize the material distribution of a cantilevered 
wing (maximization of torsional divergence stability under a mass constraint). Linear, parabolic, and piecewise 
material grading distributions along the wing span are all considered. Though this paper uses fiber volume fractions 
of composite materials (rather than metallic grading, which is the emphasis of the current work), it is notable as one 
of very few papers that consider subsonic aeroelastic tailoring via material grading. 

An FGM may be numerically realized by identifying two materials, and then specifying the spatial distribution 
of material fraction (any continuous value between 0, which is 100% of material A, and unity, which is 100% of 
material B), by which the local material properties can be tuned within the wing structure via the rule of mixtures. 
In conventional aircraft structures, stiffness can be increased using larger material thicknesses, which will also add 
mass. An FGM composed of two materials with similar density but different elastic modulus enables the 
distribution of stiffness to be tailored without affecting the mass distribution. This allows for designs with stiffness 
and mass distributions that would not be possible with conventional materials. Additionally, material and geometric 
changes in a structure have different effects on the stiffness. For example, the flexural rigidity of a plate is 
proportional to its thickness cubed, but varies linearly with elastic modulus. 

Given this, an important benchmark comparison is between the possible aeroelastic benefits of a) material 
grading between two materials that may or may not have the same density and b) those found for traditional 
thickness grading. Both studies, particularly the latter, must be formulated to maintain a constant wing weight, and 
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thus provide a more equitable comparison between the two. Given that the impact of thickness is stronger than that 
of modulus on the local stiffness, it may be found that for many situations traditional thickness grading outperforms 
FGMs, despite the weight penalty. Furthermore, many situations may be found where the best configuration is that 
which uses 100% of the stiffer material throughout the entire wing structure. This is a trivial result; this work 
desires to uncover design scenarios which can be improved with a spatial grading of material. Beyond the 
aforementioned thickness-FGM comparison, results will also be shown where both parameter types are varied 
simultaneously. A particularly interesting scenario is one where the thickness grading strategies are found to differ 
substantially between thickness-only tailoring (with a fixed material configuration) and simultaneous thickness and 
material grading. 

A simple plate-like wing, clamped at its root and immersed in subsonic flow, is used to assess the effect of 
various FGM-related parameters on the static and dynamic aeroelastic performance of the wing. The geometry of 
the plate-wing is taken from Ref. [8]: metallic material properties are used in place of the composite structures 
implemented in that work, though the mass of both plates is identical. Deviations from this straight-wing planform 
are also considered, in terms of different sweep angles and higher aspect ratios. Two performance metrics are 
considered: the flutter and divergence speeds (whichever is lower and more critical) and the maximum von Mises 
failure criteria in the wing resulting from the deformation under a steady flight condition at a specified angle of 
attack. Simple parameter sweeps are used to demonstrate changes in aeroelastic physics, though a genetic algorithm 
is used for more complex parameterizations. 

II. Aeroelastic Solution Methodology 

The modeling tools used to characterize the aeroelastic behavior of the plate-wings are described in this section. 
Though details are not provided here, these tools have been extensively verified against the aeroelastic solvers in 
NASTRAN. The wing structure is modeled using discrete Kirchhoff triangular (DKT) plate elements [9] with 
consistent mass matrices. The wing is fully clamped at the root (cantilevered). The properties of the functionally 
graded material are approximated using the rule of mixtures. The percentage of each material within an element is 
taken as the percentage at the element center. The Doublet Lattice Method (DLM) [10] is used to model the 
subsonic unsteady aerodynamics of the wing. For the cases considered here, the flow velocities are relatively small 
and therefore the Mach number is fixed at zero. 

The Finite Plate Spline (FPS) method [11] is used to transfer information between the DLM and finite element 
plate model in a work-equivalent manner, thus creating a hilly coupled aero -structural model. The premise of the 
FPS method is to create a plate finite element mesh that lies between the DLM and finite element discretizations and 
to transfer displacement and force information through this mesh by formulating constraints at locations that 
coincide with the DLM box centers and a subset of the structural nodes. DKT elements are also used to create the 
FPS mesh and the constraint equation is solved using the Lagrange multiplier method [9]. 

The p-k method is used to compute the flutter and divergence speed of the wing and obtain velocity-damping 
curves. To improve the approximation of the damping, an aerodynamic damping matrix is included in the p-k flutter 
equation [12]: 


f(D ' P 2 ' [M] - (^) • v ■ (f P ■ V ■ b ■ [Ai(i ■ *)]/*) + [K]-q- [A r (i ■ fc)]j • v = 0 (1) 

where [M] is the mass matrix, [K] is the stiffness matrix of the wing, [A r (i • /e)] and [A^i • /e)] are the real and 
imaginary parts of the aerodynamic influence coefficient (AIC) matrix, respectively, q is the dynamic pressure, p is 
the air density, v is the eigenvector written in terms of the degrees of freedom of the wing structure, V is the 
airspeed, b is a reference length (taken here as half of the root chord), i is the imaginary number and k is the reduced 
frequency: 


k = a) - b/V (2) 

where co is frequency (rad/s). The eigenvalue p is written as p = g + i • k, where the real part g is the aerodynamic 
damping imparted from the unsteady aerodynamics, and the imaginary part is the reduced frequency. 

A reduced order method is used to reduce the degrees of freedom in Eq. 1 to a small number of generalized 
degrees of freedom. The first several natural modes of the structure are used to reduce the order of the system [9] . 
For example, the reduced mass matrix is: 
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[M] = [Of • [M] • [O] 


(3) 


where [M] is the reduced mass matrix and [O] is a matrix where the i th column contains the i th eigenvector 
corresponding to the i th free vibration mode of the structure. Reduced stiffness and AIC matrices are obtained using 
the same method. The eigenvectors in [O] are normalized so that the reduced mass matrix becomes the identity 
matrix. The reduced stiffness matrix [ K ] is thus also diagonal, containing the squared natural frequencies of each 
free vibration mode. To obtain the roots of Eq. 1 using a general eigenvalue analysis, the p-k flutter equation is 
placed in standard first-order form: 


ifv\ n 17 01 

0 , El] 

jvJ'P'lo /J 

(<7 • [A r (i ■ k)] - [tf]) p-V b- [A t (i ■ k)]/k\ j 



(4) 


The highly-nonlinear dependence of the AICs on the reduced frequency complicates the p-k eigenvalue 
analysis. A valid solution to Eq. 4 is only obtained when the reduced frequency input into the AICs is exactly equal 
to the imaginary part of the eigenvalue: k = Im (p). These solutions are found using a non-iterative frequency 
sweeping technique [13]. For a set speed and dynamic pressure, k is incrementally increased and the corresponding 
eigenvalue p is monitored until a crossing with the unit-slope line (in the Im(p) — k space) is noted. The match- 
point solution for k = Im(p) is then obtained by linear interpolation. The real part of the eigenvalue indicates the 
damping, g , of the associated mode and positive damping indicates an unstable aeroelastic mode where oscillatory 
wing deformations become unbounded. An instability point is located whenever g = 0 with increased velocity. If 
the reduced frequency at this point is also equal to zero, then this point represents static divergence. Otherwise, the 
instability is a flutter point. A critical aspect of the frequency sweep approach is the tracking of eigenvalues as k is 
increased and also as speed is increased. Mode tracking is performed using an eigenvector orthogonality correlation 
method [14]. A correlation matrix is formed by computing scalar products between the current eigenvectors and the 
converged eigenvectors from the previous step. The values in the matrix are then used to re-correlate the modes. 

Static aeroelastic analysis is performed using the same tools as the dynamic flutter analysis; namely DKT 
elements to model the structure, the DLM to model the aerodynamics and the FPS method to transfer information 
between the structural and aerodynamic models. The AIC matrix is computed for zero reduced frequency and the 
following linear equation is solved to obtain the static deformed shape of the wing in steady flight: 


[K]-u = f + q- [A 0 ] ■ u (5) 

where u is the deformed shape of the wing in terms of the structural degrees of freedom, [d 0 ] is the AIC matrix 
evaluated for zero frequency and / is the force acting on the rigid wing set at a specified angle of attack. Reduced 
order modeling is not used for the static aeroelastic analysis and Eq. 5 is solved using the full order matrices. 

Three different plate-like wing designs are considered for the majority of this study: a straight wing (identical to 
the planform geometry in Ref. [8]), an aft swept wing and a forward swept wing. For each wing the semi-span is 
0.305 m and chord 0.0762 m. The sweep for the aft and forward swept wings is 20°. For the material graded designs 
a uniform thickness is used to meet the target mass of the wing. For the thickness graded designs, the thickness 
values are chosen such that the wing has the same mass as the material graded designs. In addition, higher aspect 
ratios are also considered, but the total planform area of the wing is always maintained to be equal to this aspect 
ratio 8 model. 

The steady flight condition is for an air density of 1.225 kg/m 3 and angle of attack of 5°. The airspeed used for 
the static aeroelastic analysis for the straight and aft swept wings is 15 m/s, whereas 10 m/s is used for the forward 
swept wing. The lower speed is used for the forward swept wing due to the occurrence of divergence (static 
instability) below 15 m/s. Two metrics are used to assess the performance of each wing design. The first is the 
critical aeroelastic speed: 


U C rit = min (U f , U d ) (6) 

where U crit is the critical aeroelastic speed, Uf the flutter speed and U d the divergence speed. A higher critical speed 
is desired to reduce the risk of the wing failure due to flutter or divergence. The second performance metric is based 
on a yield stress failure criterion: 
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fs — max ( Gvm! °Y)crit 


( 7 ) 


where f s is the failure metric, o VM is the von Mises stress and o Y is the yield stress (which is a spatial function 
dictated by the grading strategy). The ratio, (ovm!°y) is evaluated at the center of each finite element under the 
static aero-structural loading conditions. 


III. Parametric Studies 

A parametric study is used as the starting point for the investigation into using FGM for aeroelastic tailoring of 
a plate-like wing. There are two main goals of the parametric study. The first is to compare continuous with discrete 
material grading. The second is to compare continuous material grading, using an FGM, with continuous geometric 
grading, by varying the plate thickness. For the parametric study, the FGM is composed of two isotropic materials: 
an Aluminum alloy (Al-alloy) and Aluminum with Silicon Carbide particulates (AlSiC). Both materials have 
approximately the same density, Table 1, therefore, a change in material distribution modifies the local stiffness, but 
not the local mass (or inertia). This allows for designs with varying stiffness without changing the mass distribution. 
The density used for both materials during the study is the average of the two: 2784 kg/m 3 . 

Table 1. Material properties of A1 and AlSiC. 


Property 

Al-alloy 

AlSiC 

density (kg/m ) 

2768 

2800 

elastic modulus (GPa) 

74 

107 

Poisson’s ratio 

0.33 

0.33 

yield stress (MPa) 

276 

434 


The following discretizations are used during the parametric study. The wing structure mesh is composed of 20 
x 80 squares, where each square is divided into two DKT elements. The DLM mesh is composed of 16 x 32 boxes 
in the chord and span-wise directions, respectively. The FPS mesh is composed of 30 x 120 squares, again each 
square is divided into two DKT elements. A subset of 126 structural nodes (6 x 21) is chosen to formulate the 
constraints on the FPS mesh, in order to transfer information between the DLM and FEA meshes. 

The parametric study uses a simple linear grading strategies in the chord-wise, Figure 2a, and span-wise, Figure 
2b, directions. For the material distribution cases, a diagonal grading strategy is also considered, which has 100% of 
one material at the trailing edge root and leading edge tip and 1 00% of the other material at the other two corners, 
Figure 2c. The fraction of each material is then interpolated using standard bi-linear shape functions. In addition, a 
discrete material grading strategy is defined from the equivalent continuous grading using a step function: 


m d O) = 


1.0, 

ra c (x) > 0.5 


0.5, 

m c (x) = 0.5 

(8) 

0.0, 

m c (x) < 0.5 



where x is a center point of a finite element, m c (x) is the fraction of the material for a continuous grading and 
m d (x) is the equivalent material fraction for a discrete grading. 


A. Continuous vs. Discrete Material Grading 

The results for discrete and continuous material grading strategies are shown in Table 2, Table 3, and Table 4 
for the straight, aft swept and forward swept wings, respectively. In general, the comparison between continuous 
and discrete material grading for the same wing designs shows little difference in the performance metrics. The 
difference in the critical aeroelastic speed is less than 3%. The exception is the chord-wise grading from AlSiC to 
Al-alloy for the straight wing (Table 2), where the critical speed for the discrete grading is much lower than the 
equivalent continuous grading. This is because the flutter mode for the two designs is different and the lower flutter 
speed for the discrete grading design is due to the emergence of a hump mode. 
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Figure 2. Material grading across the flat plate: (a) chord-wise, (b) span-wise, (c) diagonal. 
Table 2. Straight wing: continuous versus discrete material grading. 


Material(s) 

A ->B (Table 1) 

Grading 

direction 

Grading 

type 

Ucrit (m/s) 

fs 

Al-alloy (uniform) 

None 

Uniform 

20.61 

0.507 

AlSiC (uniform) 

None 

24.78 

0.257 

Al-alloy AlSiC 

Span-wise 


21.97 

0.461 

AlSiC Al-alloy 

Span-wise 


23.46 

0.271 

Al-alloy -> AlSiC 

Chord-wise 

Continuous 

22.62 

0.373 

AlSiC -> Al-alloy 

Chord-wise 

22.95 

0.318 

Al-alloy -> AlSiC 

Diagonal 


22.72 

0.374 

AlSiC -> Al-alloy 

Diagonal 


22.84 

0.318 

Al-alloy -> AlSiC 

Span-wise 


21.28 

0.486 

AlSiC Al-alloy 

Span-wise 


23.85 

0.265 

Al-alloy -> AlSiC 

Chord-wise 

Discrete 

22.52 

0.398 

AlSiC Al-alloy 

Chord-wise 


17.74 

0.308 

Al-alloy -> AlSiC 

Diagonal 


22.59 

0.398 

AlSiC -> Al-alloy 

Diagonal 


22.85 

0.309 


Again, comparing the continuous and discrete designs, the maximum change between the stress failure metric is 
7%, except for the span- wise grading from Al-alloy to AlSiC for the forward swept wing (Table 4), where the 
discrete grading failure metric is 26.5% greater than the continuous value. In this case the discrete grading design 
has a lower bending stiffness compared with the continuous grading design. Thus, the discrete design deflects more 
under the aerodynamic loading. This in turn increases the load on the forward swept wing, due to its natural wash -in. 
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As both designs have a divergence speed near the flight speed, the difference in loading is magnified, resulting in the 
difference in the failure metric. For all three wings, the best material grading design has a lower critical speed than 
the best uniform material design, which is 100% of the stiffer AlSiC material in all cases. Also, with the exception 
of the aft swept wing, the lowest failure metric for the FGM wing designs is greater than that for the uniform AlSiC 
design. For the aft swept wing, the lowest stress failure metric is for the continuous span-wise grading from AlSiC to 
Al-alloy design, which is slightly lower than the failure metric for the uniform AlSiC material design. 

Table 3. Aft swept wing: continuous versus discrete material grading. 


Material(s) 

A ->B (Table 1) 

Grading 

direction 

Grading 

type 

Ucrit (m/s) 

fs 

Al-alloy (uniform) 

None 

Uniform 

18.33 

0.204 

AlSiC (uniform) 

None 

22.04 

0.149 

Al-alloy -> AlSiC 

Span-wise 


19.58 

0.209 

AlSiC -> Al-alloy 

Span-wise 


20.81 

0.146 

Al-alloy -> AlSiC 

Chord-wise 

Continuous 

20.07 

0.167 

AlSiC Al-alloy 

Chord-wise 

20.46 

0.179 

Al-alloy AlSiC 

Diagonal 


20.24 

0.165 

AlSiC -> Al-alloy 

Diagonal 


20.29 

0.181 

Al-alloy AlSiC 

Span-wise 


19.05 

0.204 

AlSiC -> Al-alloy 

Span-wise 


21.03 

0.148 

Al-alloy -> AlSiC 

Chord-wise 

Discrete 

19.96 

0.166 

AlSiC -> Al-alloy 

Chord-wise 


20.56 

0.181 

Al-alloy -> AlSiC 

Diagonal 


20.16 

0.166 

AlSiC Al-alloy 

Diagonal 


20.27 

0.181 


Table 4. Forward swept wing: continuous versus discrete material grading. 


Material(s) 

A ->B (Table 1) 

Grading 

direction 

Grading 

type 

Ucrit (m/s) 

fs 

Al-alloy (uniform) 

None 

Uniform 

11.00 

1.051 

AlSiC (uniform) 

None 

13.22 

0.275 

Al-alloy -> AlSiC 

Span-wise 


11.61 

0.707 

AlSiC Al-alloy 

Span-wise 


12.65 

0.314 

Al-alloy AlSiC 

Chord-wise 

Continuous 

11.98 

0.510 

AlSiC -> Al-alloy 

Chord-wise 

12.33 

0.392 

Al-alloy -> AlSiC 

Diagonal 


11.98 

0.518 

AlSiC Al-alloy 

Diagonal 


12.34 

0.387 

Al-alloy -> AlSiC 

Span-wise 


11.22 

0.895 

AlSiC Al-alloy 

Span-wise 


12.91 

0.292 

Al-alloy AlSiC 

Chord-wise 

Discrete 

11.90 

0.542 

AlSiC -> Al-alloy 

Chord-wise 


12.41 

0.377 

Al-alloy AlSiC 

Diagonal 


11.84 

0.553 

AlSiC -> Al-alloy 

Diagonal 


12.44 

0.374 


The results in these tables show little or no benefit of using an FGM composed of Al-alloy and AlSiC, 
compared with a uniform material wing made from AlSiC. However, if we consider the FGM as a reinforcement of 
the Al-alloy design using the AlSiC material, then some useful observations can be made. For all three wings span- 
wise material grading is the best for both performance metrics, when the root is reinforced with the stiffer material. 
This seems reasonable as the aerodynamic loading, bending moments and torsion induced shear all generally 
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increase towards the root. Thus, reinforcing this area would be expected to provide more benefit to the aeroelastic 
performance of the wing. Chord-wise grading with leading edge reinforcement is also a good strategy, with 
performance metrics often close to those obtained with span-wise grading. 


B. Thickness vs. Material Grading 

The two best continuous material grading strategies for the plate-like wing are now compared against similar 
geometric grading strategies using plate thickness variations. The two grading strategies considered are the span- 
wise and chord-wise grading with root or leading edge reinforcement, respectively. The material used for all 
thickness grading designs is the Al-alloy. For each wing, the performance metrics are computed for a range of 
thickness grading designs, characterized by the thickness of the plate at the root, or leading edge. The corresponding 
thickness at the tip, or trailing edge, is computed to maintain the total mass of the wing. The results for the straight, 
aft swept and forward swept wings are shown in Figure 3, Figure 4, and Figure 5 respectively. 
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Figure 3. Straight wing thickness grading, (a) Span-wise, (b) Chord-wise. 
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Figure 4. Aft swept wing thickness grading, (a) Span-wise, (b) Chord-wise. 

For the straight wing, the maximum critical aeroelastic speed achieved with the span-wise thickness grading is 
23.0 m/s (Figure 3a) compared with 23.5 m/s for a similar material grading using the FGM (Table 2). Flowever, all 
designs for span-wise thickness grading with a root thickness of 0.56 mm or more have a lower stress failure metric 
than the equivalent material grading design. For chord-wise thickness grading, the straight wing has better 
performance for both metrics when the leading edge thickness is 0.74 mm or more. Flowever, these designs are quite 
extreme, as the trailing edge is very thin, with a thickness of 0.14 mm or below. These results suggest that there may 
be situations where a combination of material and thickness grading would yield a better design in at least one of the 
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performance metrics, compared with a uniform thickness or uniform material design. The discontinuity in the 
critical aeroelastic speed for the chord- wise grading (Figure 3b) is caused by the emergence of a hump mode, which 
significantly reduces the flutter speed. 

For the aft swept wing with span- wise thickness grading the large drop in the critical aeroelastic speed between 
root thickness values of 0.74 mm to 0.76 mm is due to a switching of the critical flutter mode from 2 nd bending to 1 st 
torsion, Figure 4a. However, designs with a root thickness between 0.64 mm and 0.76 mm have better performance 
metrics compared with the equivalent material grading designs, Table 3. For an aft swept wing with chord- wise 
thickness grading, designs with a leading edge thickness of 0.56 mm or more have better performance metrics than 
the equivalent material grading designs, Figure 4b and Table 3. Span- wise thickness grading for the forward swept 
wing with a root thickness of 0.56 mm or more has better performance with respect to both metrics, (Figure 5a), 
compared with the equivalent material grading design using the FGM, Table 4. For the chord- wise grading, designs 
with a leading edge thickness of 0.68 mm or more have better performance compared with the equivalent material 
grading design, as seen in Figure 5b and Table 4. This suggests that thickness grading will out-perform material 
grading for a forward swept wing. 




Figure 5. Forward swept wing thickness grading, (a) Span-wise, (b) Chord-wise. 

Comparing span-wise to chord-wise thickness grading, these figures show that span-wise grading can produce 
designs with lower failure metric, but chord-wise designs can have a higher critical aeroelastic speed. There appear 
to be three main methods to reduce the failure metric. The first is to move the elastic axis (EA) forward. For the 
forward swept wing, this reduces the natural wash-in effect (leading edge up) of the wing under static aerodynamic 
loading. Therefore, the local angle of attack is reduced towards the wing tip as the EA is moved forward. This in 
turn reduces the total load on the wing. The wash-in effect is also often present for an aerodynamically loaded 
straight wing, as the EA is usually behind the aerodynamic center. For the aft swept wing, moving the EA forward 
increases the natural wash-out (leading edge down), which also helps reduce the load [15]. The load alleviation 
effect caused by moving the EA forward contributes to the reduction of the failure metric when increasing the 
leading edge thickness, or using a stiffer material at the leading edge with chord-wise grading. The failure metric 
can also be reduced using two forms of reinforcement. The first is to use a stiffer material to oppose bending (hence 
reduce the magnitude of the wash-in effect) and the second is to use a material with higher yield stress in high stress 
areas, especially near the wing root. However, reinforcement by using higher yield stress material is not necessarily 
an aeroelastic tailoring effect, in that the bend-twist coupling of the wing is not explicitly altered. 

The critical aeroelastic speed for the straight and aft swept wings is flutter, rather than a static divergence. The 
success of chord-wise thickness grading in increasing the flutter speed can be attributed to a mass balancing effect, 
as the center of gravity (CG) of the wing is moved forward [16]. Another strategy for increasing the flutter speed is 
to add mass to the tip in order to separate the natural frequencies associated with the bending and torsion modes, 
which are often involved with the flutter mode. 

The overall comparison of thickness grading with material grading for aeroelastic tailoring shows that thickness 
grading can produce designs with a lower stress failure metric or higher critical aeroelastic speed, compared with 
equivalent material grading design using a FGM composed of AlSiC and Al-alloy. In part this can be explained by 
the flexural rigidity of a plate being proportional to thickness cubed, but linearly proportional to elastic modulus. 
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Therefore, for this example, thickness variation allows more scope to move the EA. Also, the two materials in the 
FGM had the same density, therefore there is no scope for moving the CG or for using mass balancing when 
considering material grading only. However, there may be situations where a combination of thickness and material 
grading will produce a superior design, compared with a uniform material design using only thickness grading. This 
is investigated further in the next section. 

IV. Optimization with a Genetic Algorithm 

To investigate aeroelastic tailoring with both thickness and material grading, a genetic algorithm (GA) 
optimization tool is used. The initial investigation considers the straight, aft and forward swept wings and optimizes 
for both the stress failure metric and critical aeroelastic speed. The multi -objective optimization problem is: 

min (fs) and ma x{U crit ) 
such that : M < M t 

where M is the mass of the wing and M t is the target mass of the wing, which is the same as the wing mass used in 
the parametric study. The design variables are the fraction of the Al-alloy material and the plate thickness at each 
comer of the wing planform. The thickness and material distributions are defined from the corner values using 
standard bi-linear shape functions. However, given the thickness at three corners, the fourth value (tip leading edge) 
can be computed such that M = M t , thus eliminating the need for the mass constraint in Eq. 9. As such, only three 
thickness design variables need to be considered. The material fraction variables are constrained to be between 0 and 
1 and the thickness variables are constrained to be between 0.3 and 0.55 mm. However, the eliminated thickness 
variable at the leading edge tip could go below the lower limit in order to meet the mass constraint, although it could 
not become negative. 

A GA is computationally expensive. Therefore, a model with slightly coarser discretization than the model from 
the parametric study is used. The performance metrics obtained with the coarse model are within 5% of those 
obtained using the finer model. For the coarse model, the wing stmcture mesh is composed of 15 x 60 squares, 
where each square is divided into two DKT elements. The DLM mesh is composed of 12 x 24 boxes in the chord 
and span- wise directions, respectively. The FPS mesh is composed of 45 x 180 squares, again each square is divided 
into two DKT elements. A subset of 248 structural nodes (8 x 3 1) is chosen to formulate the constraints on the FPS 
mesh. The problem is solved using a non-sorting genetic algorithm (NSGA2) [17], which is a multi-objective GA 
that aims to produce a population of designs that have a good spread along the Pareto front. For problems involving 
higher aspect ratios, the number of elements in the span direction (for each mesh) is increased so as to retain square 
elements. 

A. Initial Optimization Results 

The optimization problem in Eq. 9 is solved three times for each wing. First the problem is solved using both 
thickness and material distribution as the design variables. Then the thickness variables are all set to 0.44 mm (in 
order to meet the mass constraint M = M t ) and only the material fraction values are used as design variables. 
Finally, the material variables are set to 100% of the AlSiC material and only the thickness variables are used as 
design variables. In each case the GA optimization is run with a population of 100 for 100 generations. The Pareto 
fronts for each optimization scenario are obtained and compared. The results are shown in Figure 6, Figure 7, and 
Figure 8 for the straight, aft swept and forward swept wings, respectively. 

For the straight wing, the Pareto front obtained when both thickness and material are used as design variables is 
coincident with the front when only thickness is considered. This is not surprising, as all the optimal designs on the 
Pareto front are made from 100% of the stiffer AlSiC material. The trade-off along the front is between chord- wise 
thickness grading and span-wise thickness grading. The design with the highest critical aeroelastic speed has chord- 
wise grading, whereas the design with the lowest stress failure metric has span-wise grading. This correlates well 
with the results from the parametric study in Figure 3. The optimal design when only the material distribution is 
considered is also made from 100% of the AlSiC material and the Pareto front collapses to a single design point. 
These results show that, within the scope of this study, there is no potential for aeroelastic tailoring of the straight 
wing when using a FGM composed of the Al-alloy and AlSiC materials via a bilinear distribution. 

For the aft swept wing the Pareto fronts obtained when both thickness and material are used as design variables 
is partly coincident with the front when only thickness is considered, Figure 7. However, in contrast to the straight 
wing, the Pareto front that uses both material and thickness as design variables extends further to the left and 
contains optimal designs that use the FGM. The knee in the Pareto front is due to the thickness variables hitting their 
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side limits (variable bounds). Thus, any further failure metric reduction from thickness grading is stopped and the 
optimizer found further reduction from material grading instead. The designs with the FGM are in the extended left 
part of the Pareto front and have a lower critical aeroelastic speed, by up to 3%, but a slightly lower failure metric, 
by approximately 1%, compared with the designs to the right of the knee, which are made from 100% of the AlSiC 
material. The FGM designs all have the same thickness grading, but there is an increasing amount of the Al-alloy 
material at the tip trailing edge that corresponds to the reduction in the failure metric. This reduction is due to a load 
alleviation effect, as the elastic axis near the tip is moved forward by the introduction of the less stiff A1 -Alloy 
material at the tip trailing edge. Therefore, the natural wash-out of the aft swept wing is increased, which results in a 
decrease in angle of attack near the tip and a decrease in the loading. The load alleviation effect of the natural wash- 
out for the aft swept wing can also be seen in the lower failure metric for the aft swept wing designs compared with 
the straight and forward swept wings designs. 



Figure 6. Straight plate-wing Pareto fronts. 



Figure 7. Aft swept plate-wing Pareto fronts. 

Figure 9 shows that the tip twist under the static aerodynamic loading varies linearly with the fraction of the Al- 
alloy material at the tip trailing edge, where a positive tip twist corresponds to leading edge up and increased local 
angle of attack. This demonstrates that, for an aft swept wing, a FGM composed of Al-alloy and AlSiC can be used 
as part of an aeroelastic tailoring strategy. However, the optimal designs along the Pareto front obtained when only 
material distribution is used as the design variable all have worse performance metrics compared with the optimal 
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designs when thickness is a design variable (Figure 7) and the load alleviation effect is small. This suggests that 
thickness variation has a greater potential benefit for aeroelastic tailoring than material distribution. 



Figure 8. Forward swept plate wing, all designs. 



Figure 9. Aft swept wing, tip twist under static aerodynamic loading, relating to Figure 7. 

The results for the forward swept wing are less interesting, as the Pareto front for each scenario collapses to a 
single optimal design, which is the stiffest possible design available (Figure 8) and each optimal design is composed 
of 100% of the stiffer AlSiC material. Note that Figure 8 shows the performance of all designs analyzed during the 
optimization. When the thickness distribution is a design variable, the thickness values for the optimum designs are 
maximum at the root and tip trailing edge and minimum at the tip leading edge. This result is not surprising, as the 
critical aeroelastic failure mode for the forward swept wing is divergence, which is a static mechanism and therefore 
often behaves in a similar way to the static aeroelastic solution. This is highlighted in Figure 8, which shows that in 
general, design changes that improve one performance metric, also improve the other. 

The optimization is rerun using a discrete material distribution, where the discrete distribution is defined from 
the equivalent continuous grading using Eq. 8. For the straight and forward swept wings there is no difference 
compared with the results obtained using continuous material grading. This is because the optimal designs all have 
100% of the AlSiC material and did not use any material grading. The comparison of discrete and continuous 
material grading strategies for the aft swept wing is shown in Figure 10. The discrete distribution designs do show 
some advantage in lowering the failure metric, a further 8% reduction compared with continuous grading, but there 
is no improvement in the maximum available critical aeroelastic speed through the design space. Also, for the same 
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value of the failure metric a slight improvement (~0.5%) in the critical aeroelastic speed is obtained along part of the 
Pareto front. In both cases, the improved designs take advantage of the discrete grading strategy to place the less 
stiff Al-alloy material along the trailing edge. Thus, the optimizer finds benefit from locally controlling the stiffness 
of the wing when using the discrete material grading strategy. This is in contrast to the continuous grading strategy, 
where any change in the material design variables will affect the properties at every point in the wing: an important 
difference between the grading strategies considered in this work. This suggests that the effect of choosing discrete 
or continuous grading will depend on whether local or global effects have more impact on the performance metrics. 
The realization of discrete material grading in practice can be enabled by a FGM strategy where the discrete jump is 
smoothed over a short distance, or via a mechanical weld of two different pieces. 



Figure 10. Aft swept wing, continuous versus discrete material distribution. 

B. Effect of Aspect Ratio 

The wing designs considered so far have an aspect ratio (AR) of 8. However, the use of FGM may become 
more or less beneficial when the aspect ratio is increased. The effect of wing aspect ratio is investigated for the 
straight and aft swept wings. The optimization problem (Eq. 9) is again solved using the NSGA2 optimization 
method, as discussed above. In each case both material and thickness distributions are considered as design variables 
during the optimization. If using FGM is beneficial for aeroelastic tailoring, then designs along the Pareto fronts 
should contain significant amount of both the Al-alloy and AlSiC materials. The higher aspect ratio wings are 
designed to maintain the total area of the wing. Therefore, as the aspect ratio is increased the chord length is 
decreased accordingly. The results for the straight and aft swept wings are shown in Figure 11 and Figure 12, 
respectively. 

For the straight wing with an aspect ratio of 10 or above the Pareto front collapses to a single design point, 
Figure 11. In each case the optimal design is composed of 100% of the AlSiC material and has a thickness grading 
in the span-wise direction, with a thicker root. Therefore, increasing the aspect ratio for the straight wing does not 
reveal any potential benefit of using the FGM for aeroelastic tailoring. Also, as the AR is increased, the best possible 
designs, or position of the Pareto front, becomes worse, as the wings become less stiff. 

The Pareto fronts for the aft swept wing with ARs 8, 10 and 12 contain designs that use both materials, thus the 
optimizer finds some benefit of using the FGM, Figure 12. The optimal designs for AR 8 were previously discussed, 
where there is an increasing amount of the Al-alloy material at the tip trailing edge to reduce the failure metric by 
load alleviation. This trend is also seen for the AR 10 Pareto front. However, when the percentage of Al-alloy at the 
tip trailing edge reaches 100%, the Al-alloy material is also introduced at the tip leading edge, such that the optimal 
design with the lowest failure metric has span-wise grading for both thickness and material. This design has a failure 
metric 2.2% lower than the best achieved with a uniform material. The reduction in failure metric again appears to 
be a load alleviation effect, as the natural wash-out of the aft swept wing is increased by the reduction in bending 
stiffness, resulting in a reduction in the angle of attack moving outboard. This type of load alleviation effect would 
not be seen if the loading is prescribed without considering aeroelastic feedback, and demonstrates the importance of 
including aeroelastic effects in wing structure design. 
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Figure 11. Straight plate-wing Pareto fronts with different aspect ratios. 



Figure 12. Aft swept plate-wing Pareto fronts with different aspect ratios. 

The trends observed for the AR 10 wing are also observed for the AR 12 wing. However, the introduction of the 
Al-alloy material occurs much earlier and there is only one design with 100% of the AlSiC material, which also has 
the highest critical aeroelastic speed. This is because there is no trade-off in thickness grading and all designs on the 
Pareto front have chord-wise thickness grading. The design at the left of the Pareto front has a failure metric 5.4% 
lower than the uniform material design at the right of the Pareto front. For the AR 14 wing the trend breaks down 
and the Pareto front collapses to a single optimal design that is made from 100% of the AlSiC material and has 
chord-wise thickness grading. At this point the advantage gained by allowing the wing to bend more and thus 
increase the load alleviation effect of the natural wash-out is overtaken by the disadvantage of the decreased bending 
stiffness. This shows that there is a limit to the load alleviation effect obtained by reducing the bending stiffness to 
increase wash-out. These results show that, for an aft swept wing up to a certain aspect ratio, FGM material using 
Al-alloy and AlSiC can be used as part of an aeroelastic tailoring strategy to reduce the failure metric by load 
alleviation. The results also demonstrate that aspect ratio is an important consideration when using FGM for 
aeroelastic tailoring. 

C. Different Material Combinations 

So far, this study has only considered a FGM composed of two materials with similar density, but one material 
has a 45% increase in elastic modulus and a 57% increase in yield stress, Table 1. To begin to understand the 
implications of using other materials, with more diverse properties, two additional materials are considered. These 
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are a steel material and a titanium alloy (Ti-alloy) whose properties are summarized in Table 5. These materials are 
used to create two new material pairs for the FGM. The first used the AlSiC material and the Ti-alloy. These two 
materials have similar elastic moduli (107 and 114 GPa respectively), but the Ti-alloy has 60% lower density and 
approximately double the yield stress. The second FGM material pair is Al-alloy and steel. These two materials are 
chosen as they have the greatest range in both modulus and density. 

Table 5. Additional material properties used for grading studies. 


Property 

Steel 

Ti-alloy 

density (kg/m ) 

7833 

4430 

elastic modulus (GPa) 

200 

114 

Poisson’s ratio 

0.32 

0.34 

yield stress (MPa) 

517 

880 


The density of the steel and Ti-alloy materials is significantly greater than the density of the aluminum 
materials. Therefore, the lower limit for the thickness variables is decreased to 0.1 mm. Also, due to the disparity in 
material density, it is difficult to use one of the thickness variables to eliminate the wing mass constraint in the 
optimization problem Eq. 9, when using both thickness and material as the design variables. For these cases, the 
constraint is actively used in the optimization algorithm, and all four thickness variables are considered, alongside 
the four material fraction variables. Also, when optimizing for the material distribution only, a uniform thickness 
plate is used, where the thickness is computed to meet the mass constraint. The mass constraint is the same as that 
for the initial optimization study discussed above. A series of Pareto fronts are found for the straight and aft swept 
wings at AR 8 and using different design variables. 

The Pareto fronts for the straight wing using the AlSiC / Ti-alloy FGM are shown in Figure 13. The first 
observation is that the designs that use only Ti-alloy have much worse performance than the optimal design that uses 
only AlSiC. This is because the density of the Ti-alloy is greater than that of the AlSiC material, but their stiffness 
values are similar. Thus, to meet the mass constraint the Ti-alloy designs are thinner and therefore have significantly 
lower flexural rigidity, which is proportional to thickness cubed. The Pareto front for thickness optimization using 
only AlSiC collapses to a single design point. This design has near minimum thickness at the tip trailing edge, but 
maximum thickness at the other four corners, Figure 14b. This design lies on the Pareto front obtained when both 
material and thickness are considered as design variables. However, the Pareto front has designs with slightly higher 
critical aeroelastic speeds and lower failure metric. This is due to the addition of Ti-alloy through use of the FGM. 
The design with the highest critical aeroelastic speed has a small amount (10%) of Ti-alloy at the tip leading edge 
and similar thickness grading to the optimal AlSiC design, Figure 14a. The divergence and flutter speeds for these 
designs are very close, within one speed increment (0.2 m/s) in the p-k analysis. The design with the lowest failure 
metric has an increasing amount of Ti-alloy towards the root leading edge and decreased thickness at the tip leading 
edge, compared with the optimal AlSiC design, Figure 14c. In this case, the reduction in the failure metric is mainly 
due to the reinforcement of the high stress area near the root leading edge by the material with the higher yield stress 
and is not an aeroelastic tailoring effect. 

For the aft swept wing using the AlSiC / Ti-alloy FGM (Pareto front shown in Figure 15), the optimal design 
using only AlSiC has the highest critical aeroelastic speed of all designs, and the thickness grading is similar to the 
optimal design for the straight wing, Figure 16a. Again, the optimal AlSiC design is part of the Pareto front obtained 
when both thickness and material are used as the design variables. In this case, the introduction of the Ti-alloy 
material using the FGM allows designs with a lower failure metric compared with the optimal uniform material 
design. This is achieved mainly through increasing the amount of Ti-alloy at the root, then leading edge, combined 
with a decrease in the trailing edge thickness, Figure 16b and c. The introduction of the Ti-alloy at the root helps 
reduce the failure metric as it has a higher yield stress than the AlSiC material and is not an aeroelastic tailoring 
effect. However, the accompanying thickness reduction at the trailing edge moves the CG forward, which helps 
increase the critical aeroelastic speed. 

The optimization results that used the AlSiC and Ti-alloy materials demonstrate that an FGM composed of 
these two materials can improve the performance of both straight and aft swept wings. This performance increase 
from the FGM is mainly for a reduction in the failure metric through reinforcement of the high stress areas using the 
higher yield stress material and is not an aeroelastic tailoring effect. 
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Figure 13. Straight plate-wing GA optimization using an AlSiC / Ti-alloy FGM. 



Figure 14. Straight plate-wing optimal designs using AlSiC / Ti-alloy FGM, a) best flutter design, b) best 
AlSiC material design, c) lowest failure design. Color shading shows fraction of AlSiC. 



Figure 15. Aft swept plate-wing GA optimization using an AlSiC / Ti-alloy FGM. 

The Pareto fronts for the straight wing using the Al-alloy/steel FGM are shown in Figure 17. The performance 
of the uniform material designs that used the steel material are not shown, as the optimizer fails to find any designs 
that are stable during the static aeroelastic analysis. These designs are much thinner than the Al-alloy designs, as 
steel is much denser and there is a mass constraint that needs to be satisfied. Steel does have a higher elastic 
modulus than Al-alloy, but flexural rigidity is proportional to thickness cubed and only varies linearly with modulus. 
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Thus, the steel designs were significantly less stiff and their critical aeroelastic speed is always lower than the flight 
speed of 15 m/s. The optimal design using the Al-alloy material, Figure 18b, lies near the Pareto front obtained 
when both thickness and material are considered. However, in contrast to previous results, the optimal uniform 
material design does not lie at the extreme right of the Pareto front and designs that used the FGM attained either a 
higher critical aeroelastic speed, or lower failure metric. Designs with higher critical aeroelastic speed have a 
significant amount (around 50 - 60%) of the steel material at the root leading edge, less (0-15%) at the tip leading 
edge and chord- wise thickness grading, Figure 18a. 



Figure 16. Aft swept plate-wing optimal designs using AlSiC / Ti-alloy FGM, a) best flutter design, b) design 
with U crit = 20 m/s, c) lowest failure design. Color shading shows fraction of AlSiC. 
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Figure 17. Straight plate-wing GA optimization using an Al-alloy / Steel FGM. 



Figure 18. Straight plate-wing optimal designs using Al-alloy / Steel FGM, a) best flutter design, b) best Al- 
alloy design, c) lowest failure design. Color shading shows fraction of Al-alloy. 

For these designs, the critical aeroelastic speed is flutter. The difference between the 2 nd and 3 rd natural 
frequencies (which corresponds to the bending and torsion modes involved in the flutter mode) is 68.3 rad/s for the 
design with the highest flutter speed, compared with 52.7 rad/s for the optimal design using thickness grading only. 
This suggests that the introduction of the steel when using the FGM helps increase the flutter speed by separating 
these modes. Also, the introduction of the steel material along the leading edge when using the FGM has a mass 
balancing effect that moves the CG forward, which also helps increase the flutter speed compared with the uniform 
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material design [16]. Furthermore, the critical flutter mode for the uniform material design is the 2 nd bending mode, 
and the 1 st torsion mode is a suppressed hump mode, Figure 19a. This situation is reversed for the FGM design with 
the best flutter speed, where the 1 st torsion mode is critical for flutter, Figure 19b. Although these designs are 
optimal for a deterministic problem, the presence of the suppressed hump mode is not a reliable design, as a small 
change in the design or operating conditions could have a significant adverse effect on the flutter speed. 

The designs with a lower failure metric have less steel at the root leading edge (around 20%) and mostly span- 
wise thickness grading, with some chord-wise thickness grading nearer the tip, Figure 18c. The steel at the root 
reduces the failure metric by reinforcing the high stress region with a material with a higher yield stress. In addition, 
the chord-wise thickness grading at the tip reduces the wash-in of the straight wing, resulting in load alleviation. 
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Figure 19. U-g diagrams relating to Figure 18, a) best Al-alloy design, b) best flutter design. The first bending 
(IB), second bending (2B) and first torsion (IT) modes are indicated. 

For the aft swept wing using the Al-alloy / steel FGM, the optimal design for critical aeroelastic speed is an Al- 
alloy uniform material design, Figure 20, which has minimum thickness at the tip trailing edge and maximum 
thickness at the other three corners, Figure 21a. This design lies on the Pareto front obtained when using both 
thickness and material as design variables. The designs along the front benefit from the FGM composed of Al-alloy 
and steel. The steel material is introduced at the root leading edge, Figure 21b, and the design with the lowest failure 
metric has almost 100% steel at this location, Figure 21c. This suggests that the reduction in failure metric is 
achieved by adding material with a higher yield stress to the high stress region: not an aeroelastic tailoring effect. 



Figure 20. Aft swept plate-wing GA optimization using an Al-alloy / steel FGM. 
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The optimal designs obtained using the Al-alloy and steel materials demonstrate that a FGM composed of these 
two materials can improve the performance of both straight and aft swept wings. For the straight wing the FGM 
material could be used to improve either the critical aeroelastic speed through mass balancing or the failure metric 
though reinforcement. For the aft swept wing, on the other hand, the improvement is solely due to a reduction of the 
failure metric though reinforcement. 



Figure 21. Aft swept plate-wing optimal designs using Al-alloy / steel FGM, a) best flutter design, b) design 
with U crit = 20 m/s, c) lowest failure design. Color shading shows fraction of Al-alloy. 

D. Material Grading Excluding Yield Stress Variation 

In the above examples the reduction seen in the failure metric when using the FGM often involves a load 
alleviation mechanism due to an aeroelastic tailoring effect. However, reinforcing the high stress areas near the root 
with the higher yield stress material also helps reduce the failure metric. To investigate the effect of the aeroelastic 
tailoring in isolation, a further study is conducted where the two materials in the FGM have the same yield stress, 
but either different density or elastic modulus. The baseline material used for this study is the Al-alloy. This is 
paired with a fictitious material with double the elastic modulus and another with double the density. Although the 
double density material is fictitious, its effect could be approximated (e.g., for a wind tunnel test model) by adding 
non-structural mass to the wing. Using these two new FGMs the optimization problem in Eq. 9 is solved using the 
GA considering material grading, thickness grading and both as design variables. Only the aft swept wing is 
considered, as the results so far show this to be the most promising for using FGM for aeroelastic tailoring. 

The Pareto fronts for the FGM using Al-alloy and a double density material are shown in Figure 22. The fronts 
are close when using thickness only (using the lower density material) and both thickness and material as the design 
variables. The two fronts have similar trends for the thickness variation. The best design for the critical aeroelastic 
speed has minimum thickness at the tip trailing edge and maximum thickness at the other comers. The thickness of 
the plate is reduced at the tip leading edge and root trailing edge to reduce the failure metric. Figure 23 shows a 
general correlation between the tip twist and the failure metric. This suggests that the failure metric is partly reduced 
by an aeroelastic load alleviation effect caused by an increase of the natural wash-out of the aft swept wing. This 
aeroelastic tailoring is achieved using thickness grading, as the material distribution does not affect the stiffness or 
the strength of the wing and thus does not affect the failure metric. The designs on the Pareto front obtained using 
both material and thickness as design variables do contain some of the double density material. The double density 
material is introduced along the leading edge, primarily at the tip, but also at the root, which suggests that the slight 
increase in critical aeroelastic speed is due to a mass balancing effect, as the CG is moved forward. This mass 
balancing effect helps move the Pareto front slightly to the right and results in a small gain in the critical aeroelastic 
speed, by up to 5% for the same value of the failure metric. 

The Pareto fronts for the FGM using Al-alloy and a double elastic modulus material are shown in Figure 24. 
The front for thickness variation only (using the double stiffness material) is a single point, which has minimum 
thickness at the tip trailing edge and maximum thickness at the other comers, Figure 25a. This is also the design 
with the highest critical aeroelastic speed. The front obtained when both thickness and material are considered as 
design variables contains designs with a lower failure metric by using the FGM. The front can be split into two 
parts. The first part starts from the best critical aeroelastic speed design and introduces the lower stiffness Al-alloy 
material at the root trailing edge until it reaches nearly 100%; this point corresponds to the kink in the front, Figure 
25b. Then the lower stiffness material is introduced at all four comers until the design with the lowest failure metric 
has 100% Al-alloy at the root trailing edge and tip leading edge, 86% at the tip trailing edge and 19% at the root 
leading edge, Figure 25c. The final reduction in the failure metric using the FGM, compared with the uniform 
material design, is 27.6%. The yield stress of the two materials is the same, thus the reduction in the failure metric 
must be an aeroelastic tailoring effect. Figure 26 shows a clear correlation between the tip twist and the failure 
metric. This suggests that the failure metric is reduced by a load alleviation effect caused by an increase in the 
natural wash-out of the aft swept wing. The thickness distribution for all designs on the Pareto front is the same, 
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Figure 25, therefore the reduction in the failure metric is purely an aeroelastic tailoring effect though distribution of 
stiffness using the FGM. 



Figure 22. Aft swept plate-wing Pareto fronts using an Al-alloy and a double density material. 



Figure 23. Aft swept plate-wing, tip twist under static aerodynamic loading, relating to Figure 22. 

It is also useful to compare the aft swept wing result using the Al-alloy / double modulus FGM, Figure 24, with 
that previously obtained using the Al-alloy / AlSiC FGM, Figure 7. For the Al-alloy / AlSiC FGM the reduction in 
failure metric through load alleviation is due to the introduction of the less stiff Al-alloy material at the tip trailing 
edge, resulting in an overall 1% reduction compared with the best a uniform material design. This is in contrast to 
the trend for the Al- Alloy / double modulus FGM, where the less stiff material is first introduced at the root trailing 
edge and the best design attained a 27.6% reduction in the failure metric. The key difference between the two is that 
the AlSiC material has a higher yield stress than the Al-alloy material. This suggests that the benefit of reducing the 
loading through aeroelastic tailoring is countered by the reduction in yield stress when replacing the AlSiC material 
with Al-alloy at the root trailing edge. Therefore, the optimizer favors introducing the less stiff Al-alloy material at 
the tip, rather than root, trailing edge when using the Al-alloy / AlSiC FGM. 

V. Conclusions 

This paper investigates the use of material and thickness grading for aeroelastic tailoring of a cantilevered plate- 
like wing, where the basic geometry and flow conditions are largely taken from Ref. [8]. The study explores 
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different aspect ratios and sweep angles, as well as a range of different material combinations, and a study of 
discrete versus continuous material grading. A simple bilinear grading strategy across the wing planform (no 
grading through the depth of the structure) is used for both material and thickness values. A multiobjective genetic 
algorithm is used to explore the ensuing design space, where the optimal trade-off curve (Pareto front) between two 
conflicting metrics is of interest: static aeroelastic stresses, and dynamic aeroelastic flutter. 



Figure 24. Aft swept plate-wing Pareto fronts using an Al-alloy and a double modulus material. 



Figure 25. Aft swept plate-wing optimal designs using Al-alloy / double modulus material FGM, a) best flutter 
design, b) design at kink in Figure 24, c) lowest failure design. Color shading shows fraction of lower stiffness 

Al-alloy. 



Tip twist (rad.) 

Figure 26. Aft swept plate-wing, tip twist under static aerodynamic loading, relating to Figure 24. 
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The main conclusion is that an FGM can be used to improve the performance of a plate -like wing compared 
with a uniform material design. This benefit is often seen as a reduction in the failure metric. In some cases the 
reduction is due to a reinforcement of high stress areas using a material with a higher yield stress, which is not an 
aeroelastic tailoring effect. However, for aft swept wings, the failure metric is also reduced using aeroelastic load 
alleviation by modifying the stiffness to decrease the local angle of attack of the deformed wing shape under the 
static aerodynamic load condition. There is also one case for the straight wing where an FGM design has a 
significantly higher critical aeroelastic speed compared with the optimal uniform material design. This is possible as 
there is a large difference in the densities of the two materials used to form the FGM. This permits a mass balancing 
effect to separate the frequencies of bending and torsion modes involved in the flutter mode and move the CG 
forward. Both of these effects are known to help increase the flutter speed. 

However, the introduction of FGM does not always produce an improvement for the plate-like wings, and the 
potential benefit of using a FGM for aeroelastic tailoring is dependent on several factors. The most significant factor 
appears to be the choice of material pair used for the FGM, as the trends are not consistent when different materials 
are used. The wing sweep and aspect ratio are also important factors. Depending on the sweep or aspect ratio, the 
use of FGM for aeroelastic tailoring may become more or less beneficial. Note that our study did not observe 
benefits for a forward swept wing. Also, there is a small difference when comparing discrete material grading with 
bilinear material grading, which sometimes results in an improved performance for discrete grading designs. 
However, the realization of discrete material grading in practice can be enabled by a FGM strategy where the 
discrete jump is smoothed over a short distance. This may prove to be more durable than the alternative: 
mechanically welding two different materials together where stress concentrations may arise. 

This paper has demonstrated the basic framework for material and thickness tailoring of a largely academic 
problem: a cantilevered flat plate. The longer-term application of interest is a built-up wing box structure populated 
with primary (ribs and spars) and secondary (stringers and stiffeners) reinforcement. Various aeroelastic scaling 
laws can be employed to extrapolate specific results and/or design strategies to these larger structures, potentially for 
tailoring the wing skins. A more promising strategy may be to expand the existing modeling/optimization 
framework to larger wing box structures, and replicate many of the design exercises shown here. In addition, 
thickness and material tailoring can be undertaken for rib, spar, and stringer structures as well, a situation which has 
no analogy back to the flat-plate problem (which more-closely aligns with the wing skins). The interested reader is 
referred to Ref. [18] for preliminary studies related to FGM of a built-up wing box structure. 
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